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The rare earth (R) tellurides ifeTes have a crystal structure intermediate between that of J?Te2 
and i?Te3, consisting of alternating single and double Te planes sandwiched between RTe block 
layers. We have successfully grown single crystals of Nd2Tes, Sm2Tes and Gd2Tes from a self 
flux, and describe here the first evidence for charge density wave formation in these materials. The 
superlattice patterns for all three compounds are relatively complex, consisting at room temperature 
of at least two independent wavevectors. Consideration of the electronic structure indicates that to 
a large extent these wave vectors are separately associated with sheets of the Fermi surface which 
are principally derived from the single and double Te layers. 

PACS numbers: 71.18.+y,71.45.Lr,72.15.-v,79.60.-i 



I. INTRODUCTION 

Recently the two closely related families of compounds 
i?Te2 and i?Te3 (i?=rare earth elements) have attracted 
considerable attention for their low dimensional elec- 
tronic structure and associated charge density wave 
(CDW) formation. These materials are based on sin- 
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FIG. 1: (a)Schematic diagram showing the average (unmod- 
ulated) crystal structure of 7?2Tes. 6-axis is vertical in the 
figure. Solid grey lines indicate unit cell. (b)High resolution 
TEM of Gd2Te5 along [101] direction, showing alternating 
single and double Te layers. 



gle and double Te layers respectively, separated by i?Te 
block layers. Their electronic structure is especially sim- 
ple, being determined by Te p x and p y orbitals in the 
nominally square Te planar layers. In the case of i?Te3, 
large portions of the resulting quasi 2D Fermi Surface 
(FS) can be nested by a single incommensurate (IC) wave 
vector, resulting in a sharp peak in the general suscep- 
tibility x(<l)- The material suffers a CDW distortion at 
this wavevector for all members of the rare earth series 
(R = La-Nd,Sm,Gd-Tm and Y), with T c values depend- 
ing sensitively on R due to the lanthanide contraction^ 
In contrast, the maximum in x(?) f° r the related single 
layer compounds i?Te2 (R — La, Ce) is less well defined, 
and the resulting superlattice modulation varies between 
rare earths^ The CDW gap is larger in the ditelluride 
than the tritelluride (for instance, the maximum gap in 
CeTe^ is 600meV in contrast to 400meV for CeTe^ and 
although the CDW transition has not yet been identified 
in the ditelluride, transition temperatures are anticipated 
to be somewhat higher, too). 

The title compound i?2Tes has an orthorhombic struc- 
ture (Cmcm) as illustrated in Figure [1] Note that for 
this space group setting, the long b axis is perpendicu- 
lar to the Te planes, while the shorter a and c lattice 
parameters lie in the Te planes and are almost equal in 
length^ The material is intermediate between the two 
better-known families i?Te2 and i?Te3 described above, 
consisting of alternating single and double Te layers, sep- 
arated by the same i?Te blocks (Figure [1]) . As we will 
show in this paper, the electronic structure of this mate- 
rial is reminiscent of the single and double layer variants, 
essentially comprising sheets associated with each of the 
Te layers separately. The existence of this compound 
raises the question of whether separate modulation wave 
vectors might exist on the single and double Te planes 
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TABLE I: Crystal Growth Parameters 



Crystal 


Melt Composition (at.%Te) 


Temperature Profile 


Nd 2 Te 5 


92.50% 


1050-880° C 


Sm2Te5 


90.00% 


1000-920° C 


Gd 2 Te 5 


92.00% 


1050-900° C 



separately, and if so how these wave vectors might inter- 
act or compete with each other. 

Although crystals of ii^Tes have previously been grown 
from an alkali halide flux and their average structure 
reported, to date no superlattice modulation has been 
identified for this material. In this study, we describe an 
alternative method to grow high quality single crystals 
from the binary melt, and use transmission electron mi- 
croscopy (TEM) to probe the lattice modulation. We find 
that all three compounds exhibit a modulation wavevec- 
tor oriented along the c* axis with a magnitude close 
(i?=Nd,Gd) or equal (R=Sm) to 2 /3c*, similar to that 
of the tritelluride compounds. In addition, each com- 
pound exhibits at least one further set of superlattice 
peaks oriented away from the c* axis. Calculations of 
the Lindhard susceptibility show that contributions to 
x(<t) enhancements arise from sections of the Fermi sur- 
face associated separately with the single and double Te 
planes and indicate that these different wavevectors, at 
least for R=Sm and Gd, originate from CDW formation 
in the double and single Te planes respectively. 



II. CRYSTAL GROWTH 

High quality single crystals of i?2Te5 (i?=Nd,Sm,Gd) 
were grown by slow cooling a binary melt. Inspection 
of the equilibrium binary alloy phase diagrams^ reveals 
that i?2Te5 has a much narrower exposed liquidus than 
does either i?Te2 or i?Te3, corresponding to a tempera- 
ture range of less than 50° C and a melt composition that 
varies by less than 3 at.%. Hence, for each rare earth 
it has been necessary to carefully determine the precise 
melt composition and temperature profile to achieve the 
optimal growth conditions that avoid the appearance of 
second phases. For this reason, we have focused on just 
three members of the rare earth series, Nd2Tes, Sm2Te5 
and Gd 2 Te 5 . 

Elemental starting materials of rare earth metal (Ames 
MPC, 99.50% for Sm and 99.80% for Nd and Gd) and 
tellurium (Alfa Aesar, 99.9999%) were cut and placed in 
alumina crucibles and sealed in evacuated quartz tubes. 
The ampoules were placed in a furnace and ramped to 
1050°C before slowly cooling to an end temperature (Ta- 
ble I) at which they were removed from the furnace and 
the remaining flux separated from the crystals by de- 
canting in a centrifuge. The optimal melt composition 
and temperature profile varied even for the three closely 
spaced members of the rare earth series studied here, and 




FIG. 2: High resolution TEM image of Gd 2 Te 5 looking down 
the [l01]direction. The image shows regular crystal structure 
over a macroscopic length scale (more than 25nm) without 
intercalation or stacking faults in the layering along the b- 
direction. 

are listed in Table 1. The resulting crystals were gold in 
color, forming thin, malleable and micaceous plates. 

Single crystal X-ray diffraction was used to confirm the 
phase of the crystals. 8-26 scans along the (OkO) direction 
revealed clear peaks for even k with the appropriate lat- 
tice parameter, indicating well formed single crystalline 
phase of i?2Tes. In the measurement, it was observed 
that some crystals showed weak i?Te3 peaks mixed with 
very strong -R 2 Tes peaks. The i?Te3 peaks could be re- 
duced in magnitude or even caused to totally disappear 
by removing the surface layers of the crystals using ad- 
hesive tape. This thin layer of i?Te3 forms on the surface 
of the i?2Te5 crystals during the rapid cooling, while the 
remaining melt is removed by centrifuge, and is essen- 
tially a consequence of applying this growth technique 
to a material with such a small exposed liquidus in the 
phase diagram £ 

The composition of the crystals was examined by 
electron microprobe analysis (EMPA), for Sn^Tes and 
Gd2Te5. In both cases, tellurium content was determined 
to be 72 ±1 at.%, consistent with the value anticipated 
for i? 2 Tc 5 (5/7 = 71.4%). 



III. TRANSMISSION ELECTRON 
MICROSCOPY 

Cross-sectional TEM samples were prepared using a 
"sandwich method" : a small, plate-like Gd2Tes single 
crystal (^2 x 2 x 0.1 mm) was placed in between a 
small stack of spacers cleaved from a single crystal sil- 
icon wafer that snuggly fit inside a 3mm diameter quartz 



tube. The open spaces were in-filled using Epoxy (EPO- 
TEK 353ND) to hold the stack together. A disk with 
a thickness about 0.5 mm was cut from the quartz tube 
using a diamond saw. The Gd 2 Te5 single crystal was ori- 
ented such that the &-axis lay in the plane of the disk 
and the [101] direction was close to the disk normal. The 
disk then was ground and polished to about 30~50 ^m 
thickness. A VCR dimpler was used to further thin the 
disk center area to less than 30 /im. The final thinning of 
the TEM sample was performed at room temperature us- 
ing a Fischione Model 1010 ion-miller. The milling with 
an initial setting of the Ar ion guns condition started at 
23°, 5 kV and 5 mA, at room temperature till perfora- 
tion, then followed with 3 kV and 3 mA, 15° for 15 min, 
and final with 2 kV, 3 mA, 10° for 20 min. 

A Tccnai G 2 F20 STEM (point-to-point resolution: 
0.25 nm) operated at 200 kV was employed to do 
the microstructure investigation. High resolution TEM 
(HRTEM) simulation was done using the National Center 
Electron Microscopy simulation program which employs 
the multi-slice approximation. 

Large crystalline regions were separated by residual 
flux inclusions that appears continuous along the mi- 
caceaous planes. Nonetheless, the HRTEM image of 
the Gd2Te5 single crystal shows a highly perfect crystal 
structure over a large area (~400 nm 2 , Figure [5]). Image 
matching to the simulations of the [101] image (Figure[T]) 
suggest that the isolated bright spots are columns of Te 
atoms which make up the single and double layers of Te 
planes along the c-axis direction. The elongated bright 
dumbbell features are the Gd-Te pairs in the Gd-Te block 
layers. 

Electron beam diffraction was also measured at room 
temperature using a Philips CM20 FEG-TEM operating 
at 200kV in vacuum in order to determine the ac-plane 
modulation structure in k space. Samples were carefully 
cleaved to generate thin crystal pieces with thickness less 
then 30/mi, which were mounted on a copper grid. Opti- 
mal thickness for the measurement was achieved by mak- 
ing a small hole in the middle of the crystals by ion- 
milling in vacuum for a few hours. Electron beams at 
200kV were aligned normal to the ac plane in [010] zone 
axis and selected area diffraction patterns (SADPs) from 
the fiat thin edge of the crystal hole were observed at 
room temperature in vacuum. 

All three compounds studied exhibit a complex set of 
superlattice peaks in the ac plane (Figure [3]). As previ- 
ously observed in other families of rare earth tellurides, 
the satellite peaks in the quadrant defined by the Bragg 
peaks for h + I = even translate equivalently by recipro- 
cal lattice wave vectors G = (h,k,l), h + l = even£&*&& 
The relative satellite peak positions in the first quadrant 
are listed in Table II in units of the reciprocal lattice 
parameters. 

The SADPs for Nd 2 Te 5 , Sm 2 Te 5 and Gd 2 Te 5 are all 
different, but nevertheless have some common features. 
In particular, all three compounds exhibit an 'on-axis' 
superlattice reflection <fo oriented along either the a* or 
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FIG. 3: SADPs along (010) zone axis for (a) Nd 2 Te 5 (b) 
Sni2Te5, and (c) Gd2Tes. Bragg peaks are labeled by (hi). 
Systematic modulation wavevectors are listed in Table II. 



c* direction. Within the resolution of TEM, we cannot 
distinguish a and c lattice parameters. However, high 
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TABLE II: CDW wavevectors q = aa* + f3c* . 



Crystal 


Q 


a 


13 


Nd 2 Te 5 


qo (on-axis) 
qi (off-axis) 
q2 (off-axis) 


0.000±0.003 
0.366±0.003 
0.269±0.003 


0.688±0.002 
0.269±0.003 
0.366±0.003 


Sm2Te 5 


qo (on-axis) 
qi (off- axis) 


0.000±0.004 
0.521±0.004 


0.667±0.004 
0.000±0.003 


Gd 2 Te 5 


qo (on-axis) 
qi (off-axis) 
q2 (off-axis) 


0.0000±0.003 
0.417±0.003 
0.083±0.003 


0.6871±0.002 
0.083±0.003 
0.417±0.003 



resolution X-ray diffraction on Gd 2 Te 5 indicates that q 
is, in fact, oriented along the c* direction^ For simplic- 
ity, we have also listed this lattice modulation as being 
along c* for Nd 2 Tes and Sm 2 Te5 in Table II, although 
this remains to be confirmed. This on-axis modulation 
wave vector is incommensurate for Nd 2 Te5 and Gd 2 Te 5 
with qo = 0.688 c* and 0.687 c* for the two compounds 
respectively (Figure [U[a),(c)). In contrast, the on-axis 
wave vector for Sm 2 Tes is commensurate within the res- 
olution of the measurement, with qo = 0.667c*= 2/3c* 
(Figure Hb)). 

In addition to the on-axis CDW, each of the com- 
pounds has a distinct and unique off-axis CDW struc- 
ture. 

Neglecting the small difference in a and c lattice pa- 
rameters, which is below the resolution of this measure- 
ment, the off-axis CDW in Nd 2 Tes seems to have four 
fold rotational symmetry and the lattice modulation can 
be simply characterized by a single wavevector q\. The 
second equivalent wavevector <f 2 is generated by reflection 
symmetry about the a* and c* axes. Similar symmetry 
mapping has been reported for the off-axis superlatticc 
peak q= 0.6a* + 0.2c* in LaTe^. 

The off-axis CDW in Sm 2 Tes is unique among the 
three compounds studied in that qo and q\ generate the 
off-axis higher harmonics <fi + <fo and 2<fi + <fo which are 
incommensurate in the a-direction and commensurate in 
the c-direction. All the other peaks can be expressed in 
terms of linear combinations of qo and qi as indicated 
in Figure [3]Jb) , which means that the remaining other 
peaks are higher harmonics of those two q vectors. Cor- 
respondingly, q — 2/3c* = 0.667c* and qi = 0.521a* are 
stronger in intensity than the higher harmonics q± + qo 
and 2<fi + qo. A similar higher harmonic CDW structure 
formed by linear combination of two distinct q vectors 
was also recently found in ErTe3 by high resolution X- 
ray diffraction measurements^ 

In contrast to Nd 2 Tes and Sm 2 Tc 5 , the off axis CDW 
in Gd 2 Te5 is fully commensurate in both a* and c*. 
High resolution X-ray diffraction measurements at Stan- 
ford Synchrotron Radiation Laboratory (SSRL) 10 confirm 
that all peaks can be indexed by linear combinations of 



the two wavevectors cfi=5/12 a* + 1/12 c* and <f 2 =l/12 
a* + 5/12 c*, although it is not immediately clear that 
these are really the two fundamental wavevectors since 
other linear combinations are also possible in such a fully 
commensurate satellite peak structure. 



IV. ELECTRONIC STRUCTURE 

The electronic band structure for _R 2 Tes was cal- 
culated using the linear muffin-tin orbital (LMTO) 
method within the atomic sphere approximation includ- 
ing combined-correction terms as described in RefsfTll 
and Qjl and the results are shown in Figure H] specifi- 
cally for i?=Lu (chosen to avoid the complications asso- 
ciated with the description of (band) /-electrons within 
the local density approximation). The slight difference 
in lattice parameters originating from the structural or- 




FIG. 4: (Color online) (a) Band structure of ifeTes calculated 
by LMTO for R=Lu. The blue solid lines arise from Te 5p 
states in the square planar layers, (b) and (c) The resulting 
Fermi surfaces, illustrating the trilayer splitting that arises 
due to coupling between the three Te planes. Red and yel- 
low sheets in (b) and blue and purple sheets in (c) indicate 
the weaker bilayer splitting originating from Te double layer, 
while orange and green sheets in (b) and (c) are due to Te 
single layer and the splitting is significant. The T-point at the 
centre of the zone and other high symmetry points have been 
labeled. 
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thorhombicity was ignored in the calculation and the 
value of a=c=4.36A (b/a=10.06) from Refi for Sm 2 Te 5 
was used, which is very close to a=4.34A for SmTe3.— 
Since the FS is comprised of Te bp states originating from 
the Te atoms in the square planar layers, the general 
topology of the FS is relatively insensitive to the partic- 
ular choice of rare earth atom, and indeed to changes 
in the lattice parameter of ~ 5%, allowing us to inter- 
pret these results as prototypical for all of the other rare 
earth compounds. All calculations included a basis of s, 
p, d and / states, and self-consistency was achieved at 
1280 k-points in the irreducible (l/8) th wedge of the BZ 
(corresponding to a mesh of 30 x 8 x 30 in the full BZ). 

As previously found for i?Te 2 and RTe^r^ the elec- 
tronic structure for i? 2 Te5 is two dimensional and has 
minimal dispersion perpendicular to the Te planes (Figure 
IDJb) and (c)). Six bands formed by bp orbitals from Te 
square planes were observed to cross the Fermi level (blue 
solid lines in Figure Ufa)) and the corresponding FS at 
k y — is depicted in Figure [5] 

In the tritelluride compounds, the inequivalence of the 
two Te atoms in the double square planar layer breaks the 
degeneracy of these bands, and the resulting bilayer split- 
ting has been observed directly in ARPES studies^^ For 
-R 2 Tes, in addition to this double Te sheet, there is an 
additional single Te layer, and the band structure reflects 
this via a triple splitting of its Te states. The splitting 
between states originating from the double layer is weak 
and of a similar magnitude to the bilayer splitting in the 
tritelluride compounds, whereas the splitting between ei- 
ther of these double layer states and the state due to the 
single layer is more significant 



Z N 




FIG. 5: (Color online) Fermi surface at k y = 0. Red and 
green lines approximately indicate sections of FS, formed by 
5p orbitals from the single and double Te layers respectively. 
Arrows indicate the on-axis lattice modulation q ~ 2/3c* ob- 
served in SADP for Nd 2 Te 5 , Sm 2 Te 5 and Gd 2 Te 5 . 



In addition, a small circular electron pocket due to the 
hybridization with a Lu d-state above Ep (indicated by 
red arrows in Figure |4] (a)) was observed around the T 
point. Details of this hybridization are sensitive to the 
rare earth involved in the calculation, unlike the other Te 
bp bands, and hence we might expect that the presence 
and exact volume of this pocket varies, as we progress 
through the lanthanide series. This in turn may vary the 
precise location of Ep within the Te bp bands, although 
small fluctuations of the electron pocket volume do not 
seem to significantly impact interactions between bp elec- 
trons in the Te square planes. In addition, the topology 
of this small circular section does not contribute any ap- 
preciable peak structure to the susceptibility, leading us 
to put less emphasis on it in the subsequent analysis of 
X(q). 

Even with the same a and c lattice parameters used in 
the calculation, the orthorhombicity due to the relative 
orientation of R-Te slabs in different layers is reflected in 
the electronic structure and produced unequal electron 
pockets centered at X and Z (Figure [5]). The directional 
difference in the electronic structure, in turn, suggests 
that the band splitting at the Fermi level partially de- 
pends on the relative geometry and interactions between 
Te atoms in square planes and rare earth atoms in the 
i?-Te slabs, even though the interplanar interaction is 
believed to be small. 

Nevertheless, the overall topology of the FS is surpris- 
ingly similar to the individual sections of mono- and bi- 
layer FS structures of i?Te 2 and i?Te3 (Figure 3(a) and 
Figure 5(a) in Refill!) and thus can be approximated by 
superposing the corresponding Fermi surfaces of those 
two materials with only minor changes to account for 
differences in band filling. A close investigation of the 
character of the wavefunction at each point supports this 
view and revealed that the origin of the individual FS 
sections can mostly be attributed to bp atomic orbitals 
in either the ditclluridc-like single Te planar layer or the 
tritelluride- like double Te planar layers respectively (red 
and green lines in Figure 0, if the strong orbital hy- 
bridization or orbital mixing is ignored near the band 
crossings at the Fermi level. This is somewhat as ex- 
pected, considering that the orthorhombic structure of 
i? 2 Te5 is intermediate of i?Te 2 and -RT33 with very close 
ac parameters and alternating single and double Te layers 
along the [010] direction. 



V. DISCUSSION 

A. Electron-Phonon Coupling, Lindhard 
susceptibility and CDW Formation 

CDW formation is described by the second quantized 
Frohlich Hamiltonian, 
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n,k m,k 

+ 9n,m, q al lk+q a nk (bl n _ q + b mq ) (1) 

where k and q are the electron creation opera- 
tor in the nth band, and phonon creation operator in the 
m t h mode respectively. The electron-phonon interaction 
is tuned by the coupling g n ,m,q 

between the n t h electron 
band and the m t h phonon mode and the effect on the 
lattice distortions can be shown by obtaining the renor- 
malized phonon mode and the dispersion relation from 
Equation (2): 



^ Qm.q 
Qm,q 



Pn,q 
Qm,q 



[[Qm,qi H\ , H] 



^m.q^im 



9n,m,q „ , fc Pn,q 



2w, 



V Mh 

1/2 



Mh 



Xn(o) ^ 9n } m',q 
m' 

~ UJ m,qQrn,q 

E^Qn,m,q(Jn,rn' ,q (^m,g^m' ,q) ^ •. , , 

J^Z Xn(q)Qm',q 



(2) 



where oj m , q and Q m ,q refer to the oscillation energy fre- 
quency and the Fourier component of the non- interacting 
normal coordinate of the m t ^ phonon mode respectively 
and p n ^ q indicates electron density in the n t h electron 
energy band. 

The resultant phonon mode softening strongly depends 
on the strength of gn,m,q 

and Xn(q)- While the coupling 
strength g n ,m,q singles out the electron bands and phonon 
modes relevant to lattice distortions, the distortion wave 
vectors are selected by the peak structure in \ n (q) , which 
is mainly decided by FS nesting: 



Xn(q) = Xnn'(q)+ 

"Xnn' (q) 



Xnn'% = - 



1 



(2ir) d 



^(fc+a-Aw (3) 

1BZ e n',k+q~ € n,k 



where / n (fc) and e n t refer to Fermi-Dirac function and 
the energy of the electron in n t h band. 

Although it has never been easy to obtain the exact 
m,n and q dependence of the coupling strength g theoreti- 
cally or experimentally, Xn (?) is relatively accessible from 
band structure calculations, and indeed several authors 
have argued the origin of CDW formation in both i?Te 2 
and i?Te3 in terms of simple FS nesting conditions using 
tight binding band calculations ! 2 ' 3 ' 6 ' 7 ' 8 ' 9 ' 11 ' This model 



was successful in identifying the sections of FS which 
drives the CDW modulation in these compounds, and the 
details of the nesting was found to be dependent mainly 
on the topology of the FS at the Fermi level rather than 
of the whole band structure. 

We have used the same approach for i?2Tes, calculat- 
ing the Lindhard susceptibility x(o) of the LMTO band 
structure illustrated in Figure 3] in order to examine the 
origin of the on- and off- axis CDW super lattices ob- 
served for this family of compounds. For computational 
simplicity, the two dimensional band structure at k y =0 
was considered for the summation in the 1st Brillouin 
zone(Equation 3). 

However, if we assume an isotropic coupling strength 
gn,m,q=g and consequently calculate x(o) by summing 
over all bands including inter single-double layer contri- 
butions, this quantity is found to be relatively uninfor- 
mative. Broad maxima are found centered around 0.5a* 
and 0.5c*(Figure[S{a)), but otherwise there is not a well- 
developed peak structure that would lead one to antic- 
ipate one particular wavevector to be favored over an- 
other. 

Deeper insight can be gained when more general cou- 
pling strength g's, varying for individual phonon modes 
and Te layers, are introduced. Motivated by the iden- 
tification of distinct sections of the FS associated with 
the single and double Te square planes (Figure [5]), we ac- 
cordingly divide the six bp bands crossing the Fermi level 
Ep into two relevant subgroups that form ditelluride-like 
FS sections from single Te layers, and tritelluride-like FS 
sections from double Te layers. The contributions to the 
Lindhard susceptibility from each subgroup, Xs(o) and 
Xd(q)i can be calculated as shown in Equation (4), where 
S, S' and D, D' refer to single and double layers respec- 
tively: 



Xdouble 

s 

Xsingle(q) = Xs(q) + ^ XSD(q) 
D 

XD(q) = Y XDD'(q) 

DD' 

xs(q) = Y Xss '^ 



SS' 



(4) 



Results of these calculations are shown in Figures E{b) 
and (c). Inspection of these figures indicates that Xs(o) 
and xd(q) have a more finely peaked structure than the 
total xii) (Figure Eta)). 

Using this division, and following Equation (2), the 
resulting renormalized phonon mode dispersion is given 
by: 
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(a) i 



2.9 



<D>,m,q UJ rn,q 



Mh 



XD{q) 



2.9 



<S>,m,<7 W ™:9 ^9 

; XsW ~ 



L'S m.q^ra.g 



Xcs(9) 



(5) 



This dispersion relation (Equation 5) explicitly shows 
how the phonon mode softening depends on the aver- 
age coupling strengths g<s>i 9<d> and g<sD>, and the 
Lindhard susceptibility contributions Xs{q)i Xd(q) and 
Xsd(q)i from the single, double and inter single-double 
layer contributions respectively. This division makes it 
possible, at least in principle, to identify the relative cou- 
pling strength as well as the most relevant electron bands 
from the observed lattice distortions. In the following two 
sections, we address the origin of the on- and off-axis lat- 
tice modulation with reference to these contributions. 



B. Origin of the On-axis Lattice Modulation 

All three compounds exhibit an on-axis lattice modu- 
lation with wave vector oriented along either the a* or 
c* axis. As noted in section III, TEM cannot distinguish 
these two lattice parameters, but high resolution x-ray 
diffraction experiments for Gd2Tes have determined that 
go is in fact oriented along the c* direction for this com- 
pound. The orientation of the on-axis wave vector for 
Nd2Tc5 and Sn^Tes remains to be determined, but for 
simplicity we have referred to these as also lying along the 
c* direction. The magnitude of the on-axis wave vectors 
for the three compounds are very similar, being commen- 
surate qo = 2/3 = 0.667 c* for Sm 2 Te5, with very close 
incommensurate values for Nd2Te5 and Gd2Te5 (Table 
II). 

Significantly, the on-axis wave vectors for all three 
compounds correspond to a very well-defined sharp max- 
imum in xd{<1) calculated from the sections of the FS 
associated with the double Te layers (Figure [S^b)). A 
substantial fraction of the FS coming from these double 
layers are nested by this wave vector, as indicated by 
vertical arrows in Figure [3 The actual modulation wave 
vectors are very close to the maximum in Xd(q), and 
Xd(qo) is smaller by less than 5% compared to the cal- 
culated global maximum. This striking correspondence 
is highly suggestive that the double Te layers drive the 
on-axis CDW. 

The simpler double-layer compound i?Te3 also exhibits 
an on-axis superlattice modulation, also corresponding to 
a similarly well-defined peak in the susceptibility. In that 
case, the wave vector g ^5/7c*=0.71c* over the entire 
range of the compounds (i?=La-Tm) ,i which is also very 
close to go for i?2Te5. The difference in go values between 
the two families of compounds can be attributed to differ- 
ences in band filling, as well as to the more complicated 




n 



0.85 

0.8 

0.75 



u 

n 



0.S5 

o.s 

0.75 



I 

P 



0.35 

o.a 



u 



FIG. 6: (Color online) Color scale for all panels: red high, blue 
low. (a) The Lindhard susceptibility x{<l) = Ysnn' Xnn' {<!) a t 
q y ~0, summed for all the bands including inter single-double 
layer pairs, (b) Xo(q)'- contribution to x(o) fr° m the double 
Te layers. On-axis wave vectors for all three compounds lie 
on the global maximum. Inset: a line cut following the dashed 
line, illustrating the resonant enhancement of the mixed har- 
monics go + gi and go — 2gi + 2G X for Sm 2 Te 5 . (c) xs(q)- 
contribution to x(o) from the single Te layers. Symbols rep- 
resent first harmonics of the off-axis modulation vectors for all 
three compounds, and second harmonics for Gd2Te5. Upper 
inset: a line cut following the horizontal dashed line. The off- 
axis modulation lies close to the global maximum of Xs(q) for 
Sni2Te5. Lower inset: a line cut following the vertical dashed 
line, showing the resonant enhancement of the commensurate 
off-axis lattice modulation for Gd2Tes. 
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electronic structure in it^Tes. Nevertheless, the similar- 
ity in the nesting mechanism driving the on- axis modu- 
lation in the two compounds, and its stability across the 
rare earth series for both families, is remarkable. 



C. Origin of the Off-axis Lattice Modulation 

Given the clear correlation between the on-axis mod- 
ulation wave vectors and XD(q), it is natural to reason 
that the off-axis wave vectors might be more closely as- 
sociated with the single Te planes. Indeed, with the pos- 
sible exception of Nd2Tes, none of the three compounds 
studied exhibit any obvious correlation between the off- 
axis wave vectors and Xd(q)i whereas, as we show be- 
low, there is a close correlation with Xs(o)i at least for 
Sn^Tes and Gd2Tes. However, the off-axis modulation is 
not as simple to account for as the on-axis modulation. In 
particular, there is considerable variation in the off-axis 
wave vectors for the three compounds studied (Table II) , 
reminiscent of the variation in the lattice modulation for 
the simpler single-layer compound i?Te2 (R = La - Er). 
Consequently, we consider each of the three compounds 
separately below. 

First, Sm 2 Te5. For this compound, the off-axis mod- 
ulation lies along the a-axis, with an incommensurate 
wave vector q\ = 0.521a* « 0.5a* (here we preserve the 
notation "off-axis" to indicate that q\ is not oriented par- 
allel to <fo , even though, in this case, q\ lies along a high 
symmetry direction). As can be seen from the inset to 
Figure [BJc), this wave vector lies very near to the global 
maximum in xs(q), indicating that the single Te layers 
do indeed play a significant role in driving this off-axis 
CDW. The behavior is also reminiscent of the unit cell 
doubling associated with the CDW superlattice found in 
some of the rare earth ditellurides RTe2^ 

As described in section III, Sn^Tes is unique among 
the three compounds studied in that the mixed harmon- 
ics of the CDW modulations <fo ± <Zi and <fo ± 2<fi are ev- 
ident in SADP patterns. Close inspection of the inset to 
Figure^b) reveals that these wave vectors, which have a 
different incommensurate/commensurate structure in the 
a and c directions, are in fact closely associated with no- 
ticeably significant peak structures in xd(q)- The same 
figure also shows that <fo — 2<fi + 2G(ioo), equivalent to 
<fo — 2<fi , is actually very near to another global maximum 
in xd{q) adjacent to (101). This maximum is equivalent 
to that which is close to <fo by a simple reciprocal lattice 
translation, suggestive of a resonance in the interaction 
due to the crystal symmetry: 



^9<D>,m,q=qo-2qi UJ ' m <<l='fo-'2qi 



Mh 



Xd(<Io - 2<fi) 



^9<D>,m,q=q -2q 1 LL '"l,q=qa-2q 1 



MH 



XD(qa - 2<fi + 2G(i 00 )), (6) 



XD(qo - 2(fi + 2G ( ioo)) ~ XD(qo) 



These observations suggest a significant coupling be- 
tween the two wave vectors qo and q\ in Sni2Te5. Even 
though the off-axis modulation is principally driven by 
the single Te-planes (i.e. q\ is very close to the global 
maximum in xs{q)) nevertheless, it is not insensitive to 
the double layers. In contrast, the on- and off-axis CDW 
modulations in Nd2Tes and Gd2Tes seem to be indepen- 
dent or minimally coupled to each other, without any 
commensurate / incommensurate mixing. 

In contrast to Sn^Tes, the off-axis CDW in Gd2Te5 is 
fully commensurate. Although the modulation wave vec- 
tors <fi and q2 for Gd2Te5 are different to that observed 
in Sni2Te5, both lie close to global maxima in xs{q)i sug- 
gesting that the single Te planes play the dominant role 
in driving the off-axis CDW for this compound, too. It is 
worth noting, however, that the peak structure in xs(q) 
(Figure [S](c)) is far less well developed than that in xo(q) 
(FigureJHKb)). Rather than a single global maximum with 
little in the way of additional features, xs(q) exhibits a 
broad range of maxima along sharp "ridges" (dark red 
regions in Figure E{c)) connecting four relatively sharp 
local peaks centered close to <f=(0 0.5), (0.5 0),(-0.5 
0) and (0 -0.5). The resulting figure is reminiscent 
of similar calculations for the simpler single layer com- 
pound i?Te2, which also lack a well-defined single peak, 
and for which the superlattice modulation vectors also 
vary across the rare earth series^ We will return to the 
variation in the off axis wave vectors later. 

The two commensurate CDW vectors, qi and q\, in 
Gd2Te5 span the entire commensurate CDW super lat- 
tice peaks in k space. Such an extensive commensurate 
modulation structure may not allow a simple explanation 
in terms of perturbative higher harmonics. However, it is 
interesting to note that 2q\ , a second harmonic of q\ , lies 
exactly on top of an additional weak local maximum in 
Xs(q)i which may give some hint as to the origin of this 
extensive commensurate structure (inset to Figure [Sfc)- 
the strength of this feature depends sensitively on details 
of the calculation, but appears to be robust). Specifi- 
cally, rather than just a simple perturbation of q\, 1q\ 
itself also seems to be directly coupled to the relevant 
electronic structure through the local maximum peak at 
Xs(2(fi). This effect coherently enhances the CDW insta- 
bilities at q—2qi and seems to help the commensuration 
mechanism of the off-axis CDW to extend to higher n 
harmonics or integral multiple of q\ , while it is, in con- 
trast, weakly or minimally coupled to the on-axis CDW 
leaving that incommensurate. This behavior is in dis- 
tinct contrast to that of the other two compounds stud- 
ied, neither of which exhibits higher harmonics of the 
off-axis modulation vectors. By way of comparison with 
Sm2Te5, it is also worth noting that the off-axis CDW 
peaks in Gd2Tes lie in a high but flat region, without 
any sharp peak features, when mapped onto xd{<i) of the 
ditclluride-like double Te layers. As such, and in contrast 
to the case of Sn^Tes, the off-axis CDW wave vectors ob- 
served in Gd2Te5 appear to get enhanced mainly within 
the Te single layer by this subtle interaction and develop 
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an extensive commensurate CDW structure. 

Although the off-axis lattice modulation is different for 
Sm 2 Te5 and Gd 2 Te5, and although they each exhibit dif- 
ferent resonant mechanisms which enhance the off- axis 
CDW based on interaction with the double or single Te 
planes respectively, nevertheless, the off-axis CDW for 
both compounds appears to be principally driven by the 
single Te planes. In sharp contrast, Nd 2 Te5 appears to 
defy this simple model. Specifically, the off-axis modula- 
tion wave vectors for this compound do not lie close to 
the global maximum in xs io) (triangular points in Figure 
IHlJc)). Instead, they are found near local maxima which 
have values about 20% less than the global maximum in 
both xs(q) and xd(<1), although these features are very 
weak. Taken at face value, it appears that both single 
and double planes contribute towards the off-axis CDW 
in Nd 2 Te5, although it is impossible from this analysis 
to determine whether one or the other plays a dominant 
role. 

One of the principle defining features of the off-axis lat- 
tice modulation is the huge variation between the three 
compounds studied, especially given the minimal differ- 
ences observed in the on-axis wave vector. Given that the 
electronic structure is essentially identical for all three 
compounds, this large variation suggests that differences 
in the phonon mode characteristics play an important 
role. The atomic masses of Nd, Sm and Gd differ by 
up to 10%, affecting the lattice instability through equa- 
tion 5. Presumably, the very well-defined peak feature in 
Xd{<1) (Figure [6](b)) ensures that the on-axis wave vector 
remains tied to the wave vector at which this quantity 
peaks, even as the phonon modes and electron-phonon 
coupling change. However, the more poorly defined max- 
imum in xs(<f) (Figure E^c)) is apparently not strong 
enough to completely dominate the electron-phonon cou- 
pling to the extent that variation in the phonon charac- 
teristics are able to affect the lattice modulation. This 
behavior is consistent with that of the single and dou- 
ble layer compounds i?Te 2 and i?Te3 - the former having 
a poorly defined peak in x(<z) an( i a lattice modulation 
very sensitive to changes in rare earth,— the latter having 
a very well-defined peak in x(o)i an d a lattice modulation 
that hardly changes across the entire rare earth series. 



VI. CONCLUSION 

In summary, we have presented an alternative method 
to prepare large, high-quality single crystals of i? 2 Te5 
(i?=Nd, Sm and Gd) via a binary self flux method, 
and have presented the first evidence for charge den- 
sity wave formation in this material. All three com- 
pounds exhibit an on-axis modulation with <fo m 0.68 
along the c*-direction, in combination with an off-axis 
superlattice which varies significantly between the three 
compounds studied. Based on a consideration of con- 
tributions to the Lindhard susceptibility from the single 
and double Te planes of the layered structure, it appears 
that the on-axis CDW is driven by the double Te planes, 
whereas the off-axis CDW is principally driven by the sin- 
gle Te planes. Resonant effects associated with coupling 
of higher harmonics of these modulation wave vectors 
to local features in the susceptibility of the double and 
single Te planes appear to be relevant for Sm 2 Tes and 
Gd 2 Te5 respectively, stabilizing in the first case mixed 
commensurate/incommensurate off-axis harmonics, and 
in the second an extensive commensurate structure asso- 
ciated with just the single Te layers, decoupled from the 
on-axis modulation. 
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